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We derive optimal schemes for preparation and estimation of relational degrees of 
freedom between two quantum systems. We specifically analyze the case of rotation 



(N 
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. parameters representing relative angles between elements of the SU(2) symmetry 

group. Our estimation procedure does not assume prior knowledge of the absolute 

in 

spatial orientation of the systems and as such does not require information on the 



underlying classical reference frame in which the states are prepared. 



I. INTRODUCTION 



Estimating the state of a given quantum system is a fundamental primitive of many 



quantum information tasks. This problem is usually translated to the estimation of the 



value of a physica 
many instances 



)arameter describing specific properties of the preparation procedure. In 



global parameters of the state space define a natural 
scheme for encoding quantum information. The global parameters describe collective degrees 
of freedom of a system with respect to the external environment and are often related to an 
overall symmetry transformation of the state. 

However, encoding information into global degrees of freedom may be often problematic, 
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due to lack of knowledge of the reference frame with respect to which they were prepared, or 
due to collective decoherence by which they are affected. Encoding information into relative 
degrees of freedom, possible whenever a quantum system is decomposable into parts, can 
overcome many of the difficulties encountered in these situations. Such an encoding scheme 
las been demonstrated experimentally U] and can be applied to quantum computation 



12 1, communication 



[13I . 14 1 and cryptography Q • 



The aim of this work is to develop efficient preparations and measurement schemes for 
the relative parameters describing symmetries between different components of a system. 
We note that such measurements can induce relative relations when previously absent, as in 
the case of the relative pha se between two Fock states or the relative position between two 
momentum eigenstates [171 ]. 

In this paper we specifically confront the task of efficient estimation of relative rotation 
angles between two representation vectors of the £77(2) symmetry group. This problem was 
first addressed by Bartlett et. al. [18], who explicitly worked out the estimation of £77(2) 
rotation angles between two spin coherent states. In this article we proceed along the lines 
of earlier work that Asher Peres started together with us, and propose an extension of the 
previous methods. Our approach is based on an optimization of the quantum states used 
in such protocols with respect to some average measure of success of the estimation task, 
which we shall refer to as the fidelity. The key to this problem lies in the decomposition of 
both the signal and the measurement elements in irreducible components, invariant under 
global rotation transformations. 

In the following section we discuss the general mathematical structure of the problem. 
In section IIII Al we derive the optimal measurement for the case in which one system is 
comprised of two spin-1/2, and the other of one spin-1/2. We find that preparing two spin- 
1/2 parallel to each other leads to a marginally higher fidelity than the antiparallel case; then 
we determine an optimal preparation procedure which gives a higher fidelity then the ones 
achieved by the above preparations. This is in contrast to the known results for transmitting 
a spatial direction We then proceed in section UlI Bl bv replacing the single spin state 
of the second system with a spin-? coherent state and determine the optimal preparations 
for the cases of anti-parallel and parallel spins, which now yield nearly the same fidelity 
for any value of j. We then study the quantum/classical correspondence by considering 
the limit in which j becomes very large. In section IIII CI it is shown that in this limit the 
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problem reduces to estimating the first state with respect to a classical reference direction. 
Our results establish the correspondence between relative degrees of freedom in quantized 
systems and collective degrees of freedom defined with respect to a classical reference frame. 

II. ESTIMATION OF RELATIVE PARAMETERS 

A. Formulation of the general problem and basic notations 

We begin this section by formulating the problem for a general symmetry group and 
introducing the basic notation that is useful for our general scheme. In the following, let 
G denote a symmetry group, compact or finite, that describes the global properties of the 
system through its action on a set of parameters T. We shall consider G = SU (2) acting on 
quantum spin states, parametrized by the set of rotation angles. The possible states of the 
system are pure states ^(O), O G T in a ci-dimensional Hilbert space H carrying a unitary 
representation {U(g),g G G} of G. To introduce relative symmetry transformations, we 
assume that the representation space H is a tensor product of two components H = Hi ®H 2 
of dimensions di,d 2 respectively. The representation of G on H is decomposed into the 
product {Ui(g) ® U 2 (g)} on Hi <E> H 2 . This product representation of G x G, in which 
each component is transformed by the same element g G G, is isomorphic to G itself. We 
introduce another set, t, of parameters 9, to describe a relative symmetry between the two 
components, represented by a group of transformations G, which can be the same G as 
before or a subgroup of the latter. G refers to a symmetry property of one of the subsystems 
(say 2) with respect to the other (say 1), such as, in our case, a relative rotation angle. We 
call U 2 (h) its representation operators on the space H 2 . We shall restrict ourselves to the 
case where each of the subsystems is prepared in a pure state, although the formalism can 
easily be extended to mixed states. The total state on H can be written in terms of the two 
sets of parameters as \l/(0, 9). Its transformation under a global operation is 

u(g)*(e,e) = um <g> u 2 (g)*(e,e) = *( g e,e). (i) 

The objective of the construction is to define an efficient estimation procedure for the 
relative parameters 9, overlooking the information carried by the global parameters 6. Note 
that we might be interested only in estimating a subset of the relative parameters (which will 
be the case in the following sections). In order to quantify the efficiency of our estimation 
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procedure, we choose a utility function /(/x, 9) which measures the deviation of the estimated 
parameters \i from their true values 9. We consider only utility functions which are invariant 
under global rotations. The measurement apparatus, represented by the POVM {E^}, 
should be constructed such that it maximizes the average fidelity, denoted by F, and given 
by 

F {E^} = W d9dQP(Q, 9)Tr [p(0, 9)E„] f(ji, 9), (2) 

where P(Q, 9) is a prior probability distribution over the global and relative parameters, and 
p(0, 9) = | ^(O, 8))($>(0, 9) |. As noted in 3], we can assume that the global and relative 
parameters are independent random variables, and that the global parameter is uniformly 
distributed on its domain of definition, i.e., 

P(Q,6)d9dQ = p{6)d0dQ. (3) 

Now, using the definition for the global transformation (JJJ and the properties of the trace, 



we can write 



where 



F {EJ = £ / dffp(0)Tr [p{9)E,] f(fi, 9), (4) 



p(9) = J dgU\g) P {Q,6)U{g). (5) 

In the last equation, dg is the invariant measure for the group G. As a consequence of 
equations (jlj) and (0), we need only to consider a reduced form of the input state which is 
manifestly invariant under global transformations. Schur's lemma [2^ then assures that the 
input state is block diagonal in the irreducible representations of G, 

P{9) = Y J P (J \0)- (6) 

The above considerations also have implications on the form of the optimal measurement. 
In fact, using the global invariance of p, we have 



F{E,} = F{u{g)E,U\g)} 

= f{J dgU(g)E lx U^g)Y 



This relation implies that when searching for the optimal estimation procedure it suffices to 
consider POVMs which are invariant under a global transformation 

U(g)E tl tf(g)=E tl ,\/geG. (7) 
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Equation (JJJ) is by itself a strong prerequisite on the structure of the POVM elements. Indeed, 
if combined with Schur's lemma, it implies that whenever the total Hilbert space 7i can be 
decomposed into a direct sum of irreducible representations under the global transformation, 
the optimal POVM elements labelling the different outcomes have the simple form 

E, = Y,P^ E t> ( 8 ) 
j 

where each of the operators E J has support only on the representation labelled by J and 
Y^ilP^jE^ = lj. Here the symbol lj denotes the identity operator in the J subspace. The 
search for optimal POVMs may be further restricted to the case in which all elements have 
support on only one representation space, due to the linearity of the fidelity functional (given 
a POVM of the form (jHJ), the POVM with elements E^j = p^^jE^ yields the same fidelity 
as the original one by linearity of the trace and obeys the restriction that each element has 
support on only one representation). 

In the following sections, the general considerations outlined above will be applied to the 
problem of transmitting relative rotation angles of the SU(2) group. 



B. Estimation of a relative angle between spin coherent states 

In estimating relative, as opposed to absolute, rotation angles, we assume no prior knowl- 
edge of the overall orientation of the classical frame in which the system is defined. Following 
the notation introduced above, the problem may be illustrated as follows. 

Imagine that, with no prior knowledge on the absolute spatial orientation of two observers, 
Alice and Bob, we were requested to estimate the angle /3 between two unit vectors hi and 
n 2 , each chosen by one of them, by measuring a pair of SU(2) spin states prepared in the 
corresponding reference frame of each observer. This task is in general possible owing to the 
fact that states belongingto an 577(2) representation space can be used as intrinsic direction 
indicate, |M flt I Q fl Q and therefore it mak es sense to consider rdative angle 
between them. 

The simplest way to achieve the task would be to consider two SU(2) coherent states, 
corresponding, say, to spins ji, j 2 , 



"1 • Jl|V>l) = Ml), "2 • 3 2^2) = j2\lp2)- 



(9) 
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Without loss of generality, we can assume that Alice and Bob choose the z-axis of their 
reference frame. A state denoted by 1 1P2) in Bob's frame is written in Alice's reference frame 
as C/f (a)U 2 y (fyU^ (7) l^), where U%(a) = e iJ,i ' xa , and similarly for the other directions. 
The angles a, /3,7 are the three Euler angles relating Alice's reference frame to the one of 
Bob. Note that the angle (3 G [0, n] is also the angle between fii and n 2 . We introduce a 
global reference frame which is rotated with respect to Alice's frame by an angle a around 
the z-axis. In this frame the composite state is given by 

I 7) > = Uf{a) I Vi) ® f/ 2 2 (/3)^(7) 1^2). (10) 

For spin coherent states, as in Eq. Q, the above equation is simplified to 

|*(/3)) = l^m! = j x ) <g> U 2 y ({3)\j 2 ,m 2 =j 2 ), (11) 

up to an overall phase. Notice that in Eq. (fTUj) . (fTTj) we have implicitly specified the global 
parameter B, which we shall omit from now on from our notation. So far we have a set of 
three relative parameters, 9 = {a,/?, 7}, with a joint probability distribution given by the 
Haar-measure 

p(a,0,j) = sin (3, (12) 

which corresponds to a random orientation of Alice's and Bob's reference frames. Since the 
party making the measurement (Bob) is not interested in estimating a, 7 € [0, 2tt], these 
parameters are averaged out by integrating over their range. We are then left with the 
probability distribution p(f3) = sin/5/2, which corresponds to the probability density of the 
angle between two random unit vectors in three dimensions. 

Let us denote by Hi,H.2 the Hilbert spaces of the systems prepared by Alice and Bob 
respectively, carrying the SU(2) representations j\ and j 2 . The composite Hilbert space 
H = Hi®H,2 carries a diagonal product representation of SU(2)<g>SU(2), W 1 (g)<S>U : ' 2 (g), g G 
SU(2), which corresponds to a global symmetry operation labelled by the parameter g, acting 
identically on the two subspaces. This representation may be reduced as Yj^jS^-m 7~Lji 
where each component has multiplicity one. Invariance of the measurement operators under 
a global rotation (as explained in the previous section) reduces the signal state to the form 



J 



(13) 
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where ITj are projectors on the representations J. Since each representation of the global 
rotation (specified above by J) appears only once in the signal state, the measurement 
process amounts to estimating a probability distribution over the relative angle (3. The 
scenario described above, where each of the parties prepares a spin coherent state, is the 
one examined in [181. 



III. GENERAL SCHEME 



In the procedure outlined in the last section, the two parties, Alice and Bob, use spin 
coherent states in order to indicate their chosen direction. However, for the task considered 
here this is not the optimal preparation. It is known, in fact, that optimal 577(2) direc- 
tion indicators exploit entanglement between components belonging to different irreducible 
representations ja, 0, 0, 0, Uj. This suggests to consider the following general encoding 
procedure. Let 

Jmax Jl 

l $ >= E E <\jirm) (14) 

31=0 m\=-ji 

be a generic state in Ti\. By choosing a unit vector hi, Alice would prepare the state 

|Vi> = £/(6i)|<&), (15) 

where U(hi) is a unitary operator corresponding to the rotation which carries Alice's z-axis 
onto hi. Our goal is to find the optimal state | $) for the case in which Bob indicates his 
direction h 2 with a coherent state | ip 2 ) satisfying 

h 2 ■ J 2 | -02> = 32\ tp2)- (16) 

Note that equations (JT5j) and ()16j) are written in Alice's and Bob's reference frames, respec- 
tively. In a global reference frame, specified as in the last section, the total state is given 
by 

|*(a,/3 )7 )> = Ut{a)\^) ®U y 2 {P)U z 2 {i)\^ 2 ) 

= E <e imia \3i™i)®U2(P)e ii2,r \hh)- (17) 

jl m\ 

As before, the state | is expressed in terms of an arbitrary orientation of the global 
reference frame, so that the parameter may be omitted. However, one can see that the 
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angles a, 7 now induce relative phases between the different components of the state. Since 
our protocol does not deal with the estimation of a and 7 (only the angle j3 is considered), 
we will average over them in the expression of the fidelity 

F = E / daPjdgTr [u(g)p(a, P,j)U\g)E tl ] f{^(3), (18) 

where we denoted <i a/ 3 7 = 1/87T 2 sin(/9) dadftd^. Integrating over a and 7, we have 

j dad 1P {a,(3, 1 ) =^c miPmi ®U!(f3)\j2j2)(j2j2\U! j (f3), (19) 



/Hi 



4tt 2 

with 

Pmi = IVvKVv I, |^mi)= E a^Jji,mi), (20) 



J ma 2 



and c mi given by 



Ji=mi 



Cmi=El<| 2 , E < W = 1 - ( 21 ) 
ii m i 



In the following we will search for the optimal generic state | $). Note that after integrating 
over a, 7 we get a convex combination of states with different mi. Thus the fidelity will 
contain a linear combination of contributions from the different mi sectors, 



F = ^c mi F( mi ) < maxi?(mi). (22) 

mi 

From the above discussion it is clear ^| that the optimal generic state can be taken with 
mi fixed, i. e, of the form 

\*™)= E <bi^i)- ( 23 ) 

ii=mi 

In the following we shall restrict ourselves to generic states of this form. The rotation in 
Eq. (JT7j) by the relative angle (3 6 [0, tt] is expressed in the standard Euler angle notation as 
a rotation around the |/-axis, given by the matrix 

U^)\j 2 j 2 )=J2d%M\hm), (24) 



where the d™, m ([3) can be expressed using Jacobi polynomials (see for example j2l|). The su- 
perscript in the above equation refers to the (2j 2 + l)-dimensional irreducible representation 
of spin j 2 . The signal state (fT7|) is given explicitly (up to an overall phase) by 

I *(/*)> = E<E4 2i2 (/9)b'iWi, J 2 m 2 ), (25) 
ii m 2 
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where lii m-u jp m 2 ) = \jiTni) <S> \ j2 m 2)- Note that unlike the example discussed in section 
III Bl and in [18|, we now exploit repeated irreducible representations J, since each value of 
ji gives rise to a series of total angular momentum J = \j\ — (ji +32)- The equivalent 

repetitions of the representation J are labelled by ji, with j 2 being fixed. The state can be 
written in the basis (J M,j\ j 2 ) using the unitary transformation 

\j imi , j2 m 2 ) = £q^u 2 »J JM > J1J2), (26) 

J M 

with C denoting the Clebsch-Gordan coefficients, and M denoting the z component of the 
total angular momentum. 

We now need to compute the averaged state p(/3), following Eq. A corollary to 

Schur's lemma 2Q| then states that, for irreducible representations a, r of a group G, a 
group-averaged operator satisfies 

dgUM(g)AUV(g)i= 6 ^ A)1 , (27) 

where 5 a>T is a Kronecker delta over the inequivalent representations and the trace is com- 
puted over the d a - dimensional space of the irreducible representation a. The above corollary 
is applied to obtain the invariant reduced density operator p(f3). Remembering that the 
global rotation operator U(Cl) is a direct sum of operators U(Q) = (BjU^ J \ri), we see from 
Eq. (j2*7|) that p{(3) is also block diagonal in the representations J, 

p(p) = f u{Q)p{(3)u\n)dn 



= EP (J) (/?), (28) 
where the operators p^ J \P) bear the indices and are given by 

P (J) (/?kii = 'E<J,M,j[\ptf)\J,M,j 1 ). (29) 

M 

Note that by averaging over the global rotation one does not diagonalize the operator p(/3) 
with respect to the additional quantum number j\. The invariant signal state p{f3) is block 
diagonal over the irreducible representations with off diagonal elements across the repeated 
ones 



1 2 

P{0) = X! X! ToTXTT ^ a ™i fl ™' (^?2 m2(£0) ^jim 1 j2m2^j' 1 m 1 j2rn 2 \^ h)(Ji (30) 



Remember that mi is fixed. 
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Note that the state p(/3) does not, even implicitly, depend on the orientation of the 
reference frame of the measurement apparatus. We will now use p(/3) to determine the state 
| ipi) which will enable optimal estimation of (3. To this end, we fix a convenient figure of 
merit as measure of the discrepancy between the estimated and the given value of /3 £ [0, it], 
namely the quadratic utility function /(/x, 0) = cos 2 ((/i — /3)/2), where /i is the estimated 
value of the parameter. The choice of the utility function is not unique, and a different choice 
might lead to different optimal states and POVMs. However, the optimization procedure, as 



described below, is independent of this choice. The above uti 



ity function 



nas 



;he advantage 



of having been broadly used throughout earlier literature 

We denote the POVM elements by {E^}, with p, G [0, 7r], and Y^^E^ = 1. The average 
fidelity with respect to the given figure of merit, integrating over all possible transmitted 
angles f3 and all possible inferred values p, is 

F [{p}, {E,}} = £ / Tr[p(/3)^] cos 2 (^^j sin/3 d/3/2. (31) 

Note that the fidelity is a functional both of the set of estimates {p} and of the POVM used 
for the estimation procedure {E^}, where to each estimate corresponds a 

(single) POVM element. The probability of estimating p for a true angle f3 is Tr\p(f3)E fl ]. 

The above expression can be rewritten by exchanging the order of the integral with the 
trace (due to the linearity of the integration and finiteness of the sum) and gives the fidelity 
in the form 

F [{p}, {E,}]=J2Tr{A,E,}. (32) 
n 

with 

= J p((3) cos 2 ((/i - p)/2) sin (3 dfi/2. (33) 

Since p{(3) is block diagonal, also can be written as a direct sum A^ = A^. 

In terms of the basis representation states of angular momenta (</;ji,j2), and using 
Eq. ([30)1 . the operator A^ can be explicitly written as 

Cjimi j2m 2 ^J>i j 2 m 2 } I J) 3l){J-i 3\ I ) (34) 



where 



I% 2 (P) = I K m2 ((3)) 2 cos 2 ((p-(3)/2)smf3df3/2. (35) 
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This expression may be evaluated using the properties of the Wigner functions d 3 , and 



their representation in terms of Jacobi polynomials 



2l| . The average fidelity can now be 



written as a sum of contributions from each subspace of given J 

where Y^^E^ = tj. Given a set of estimates {/i}, the task of maximizing the expression 
J2n Tr^^E 1 ^) is straightforward, at least numerically (for example, by using semidefinite 
programming j^])- However, in our approach, in order to maximize the average fidelity (|36jh 
we need to maximize over the set {/i}, so that the maximal fidelity -F max will actually be 
given by 

F max = max max F [{//}, {E^}] . (37) 
M {e^} 

A. The case of j 2 = 1/2, j 1 € {0, 1} 

We will now solve the optimization problem of Eq. (|37j) for the following case: the state 
| ip2) is a spin 1/2 coherent state, while | is composed of two spin 1/2 systems, so that 
ji G {0, 1}. According to the discussion in the previous section, we can restrict ourselves to 
two classes of generic states | $) 



| $ ) = a \j! =0m 1 = 0) + VI - a 2 \j x = lm 1 = 0), (38) 

and 

|$i) = \ji = lm 1 = 1) (39) 

Let us first discuss the case where | <3>) = | In this simple case, the state | ipi) is just 
a spin-1 coherent state or, viewed as composed of two spins, it is a polarized state with 
parallel spins along the vector ni. Coupling the representations of ipi and ip2 gives 

1® 1/2 = 1/2 3/2, 

and the operators are one-dimensional, making the optimization trivial. In this case, 
the optimal measurement simply consists in the projections onto the J = 1/2 and J = 3/2 
subspaces, as there are no repeated representations. The fidelity achieved with this state is 
F max [parallel] = 0.90983. 
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Next, we consider | $) = | $ )- Coupling the representations of ip\ and ip2 gives in this 
case 

(1©0)® 1/2 = 1/2 © 1/2 © 3/2, 

and therefore the density matrix p{(3) contains two blocks of dimensions 2 and 1 correspond- 
ing to J = 1/2 and J = 3/2, respectively. The operators are given by 



( a 2 (4+7rsin^t) ay/1— a 2 cos fi 

8 6^ 

ay / I^"cos At 6(2-a 2 )+37r(l-a 2 )sm M 

6^ 72 

and 

A 3/2 = 12 (1 - a 2 ) + 3(1 -a 2 ) sin // 
„ 36 ' 1 J 

We are seeking the set /i and which maximizes the mean fidelity (jHEJ)- Let us start with 

the J = 3/2 subspace. Since this subspace is one-dimensional, the restriction to operators 

which are invariant under a global rotation leaves us with one operator only, E^ , which is 

the projection operator on the J = 3/2 subspace. The estimate /i 3 / 2 which maximizes the 

corresponding expression for A 3 J 2 in reference to Eq. ([40)1 . is obviously given by ^3/2 = tt/2. 

Next, we consider the 2-dimensional subspace of J = 1/2. Following j^j], we define an 

operator T as 

t=j:v^ ( 4i ) 

For a set {/i}, a POVM {E^} is optimal if and only if it satisfies the following set of conditions 

T - > (42) 

for each /1 in the set of estimates {^t}, with the additional requirement that T be hermitian. 
The inequality sign in Eq. (|4*2*|) means that the operator T — A^ must be positive semi- 
definite. The maximal fidelity will then be given by 

F max = TrT. (43) 

In order to see that equations ()42)1 and ()41)1 indeed lead to the maximization of the mean 
fidelity, consider a different POVM {E^}, such that 

EK = 1 - ( 44 ) 
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The difference between the fidelity achieved with this POVM and the one achieved with the 
optimal one is 

F max - F' = Tr ]T(T - A,)E'^ (45) 

thanks to Eqs. (|4*H|) and (jUJ). Now, if C and .D are positive semi-definite hermitian operators, 
then they satisfy 

Tr(OD) > 0. (46) 
Setting C = T — and D = E'^, we obtain 

F max - F' > 0, (47) 

as desired. 

Let us first maximize the average fidelity for only two estimates /ii and /i2, in correspon- 
dence to which we have the POVM elements Ej/ 2 + E 1 ^ 2 = lj=i/2- Then 

= « - > 0, (48) 

since T — A 1 ^ is non-negative if the POVM is optimal. Let us denote by rji and |^) the 
eigenvalues and corresponding eigenvectors of the operator A = A 1 ^ — A 1 ^ 2 . For each \r]i) 
we can write 

(rkU% - A%)E%\ Vi ) = th(th\E£\th) > 0, (49) 
using (}4*%j) . If we assume that rji is negative, then (J4T?)) gives 

fal^l**) < 0; 

on the other hand, since Ej/ 2 is positive semi-definite, we must have 

(rulEgh) = 0, if Vi < 0, 

and similarly 

(VilE^lvi) = 0, if Tfc >0. 

Thus i?^ 2 projects onto the subspace spanned by \r]i) with rji < and .E^ 2 projects onto 
the subspace of positive eigenvalues of the operator A. The subspace with 77^ = does not 
contribute to the fidelity, so that the maximal fidelity is given by 

F max = Tr T = Tr A 1 '* + £ Vi . (50) 

r H >0 
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Let us now assume that fi\ = // and /X2 = tt — //. We shall see that this choice will lead to 
the optimal measurement for the class of states under consideration. Indeed, we have 

' g a Vl— a 2 cos 11 



A 



3^3 

a\/l—a' 2 cos fi ri 
V 373 U 



(51) 



with eigenvalues ^v* 1 ^<^J 2EJ± an( j corresponding eigenvectors |+) = (1, 1) T and |— ) = 
(1, — 1) T . The contribution to the fidelity from the J = 1/2 subspace F 1/2 is now given 
by 

F"'M = aVT ^ CO! " -^A>>°. (52) 

At this point, in order to find the maximal mean fidelity under our assumptions, it suffices 
to maximize the function F 1/2 (fi). A simple calculation shows that the maximum is attained 
for 

v = tan" 1 [3^(1 + 2a 2 )vr / (8aVl - a 2 )] . (53) 
It remains to check that indeed the choice 

V\ = V, /jL 2 — 7T - V (54) 

leads to the maximal fidelity for the J = 1/2 subspace. A proof of this fact is provided by 
the following argument. Consider a general set of estimates {/x} = /x 2 , 1^3, ...,fi n }, with 
a corresponding set of POVM elements Ej/ 2 , and let 

^[{ M }]=maxF^[M,{^}] (55) 

be the maximal fidelity achieved by this set. We would like to show that by adding v and 
71 — v to the set {/x} the mean fidelity can never decrease with respect to the optimal bound 
and, at the same time, the bound is attained by these two values alone. Let {ft} denote 
the new set obtained by adding v and n — u, defined by Eq. to the set {/it}. The first 
property, 

F 1/2 [{£}] > F^ 2 [M] , (56) 

simply follows from the fact that adding estimates to a given set can only increase the mean 
fidelity. To complete the argument we still have to show that 

F^ 2 m] = F^ 2 [{u,n-u}]. (57) 
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The optimal measurement for {u, tc — z/}, following the earlier discussion, is defined by the 
projectors on the negative and positive eigenvectors of the operator A]/ 2 — A^'l^, i.e., 

K /2 = \-)H, E£ v =\+)(+\. (58) 

Consider a POVM for the set {/2} which consists of the two operators in ([58)1 . and of 

Ej/ 2 = if il^v^-v. 

This POVM is optimal also for the set {/i}. To see this, we need to check whether the 
condition 

T - A 1 ' 2 > 

holds for all // € {/2}, with 

T = Ay\-)(-\+A£ v \+)(+\. (59) 

Let us evaluate the entries of the operators T — A 1 ^ 2 in the basis |+), |— ). These are given 
by 

(-|T-AV 2 |-) = (-K /2 |-)-(-|^ /2 |-) 
(+|T-^ 2 |+) = (-|^ 2 |-)-(+l^ /2 |+) 

(+|T - Air\-) = (+\AT\-) - (+\A^ 2 \-), (60) 

where we have used (+|Aj r / f i ,|+) = ( — \Al /2 \—). The eigenvalues Ai(/i), A 2 (/i) of the operator 
T — A 1 ^ 2 can now be calculated from ([60|h and their positivity can be verified (at least 
numerically). The positivity of these eigenvalues (of which we do not report here the explicit 
expression) implies then (|H7j) . For all input states \ipi) discussed in this paper we have 
verified that the POVM given in (|58)1 . with v given by (|5T?|) . is indeed optimal. 

From the above optimization procedure we see that the maximal fidelity, as a function of 
the parameter a, is given by 

F_ K = 0] = aVT ^ C ° SV + Tr AT + T. A%1, (61) 

with v given by (fo^j) . The fidelity as a function of the state parameter a is plotted in 
Fig.HJ 
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FIG. 1: Fidelity as a function of the state parameter a. 

The maximal fidelity is achieved in correspondence of the state | ip op t), by setting a = 0.609, 
and is F max [i(: opt \ = 0.91092. For comparison, the anti-parallel spin state 

Isanti) = | U) = i|00) + i|10) 

leads to a fidelity of .F max [| ||)] = 0.90982, lower than the one obtained using the parallel 
spin state by only a factor ~ 10~ 5 . Similar results for parallel and anti-parallel spin states 
were obtained by N. Gisin and S. Iblisdir |24j . 

To conclude this part, we compare the above results to earlier results on quantum direction 
indicators, where a quantum system carrying a representation of the rotation group is used 
to transmit a spatial direction between two observes that do not share a common reference 
f ram e. A S show, in QQ, if the state of the q nnntn m syS te m is constructed f rom two spin 
1/2, (i.e., constrained to have maximal spin 1), encoding the directional information into 
anti-parallel spins proves to be the optimal strategy. Here we see instead that if the receiver 
is interested only in the relative orientation of this state with respect to another state, the 
anti-parallel spin state gives nearly the same fidelity as the parallel one, which is well below 
the optimum. 



B. Higher values of j2 

The estimation of the relative orientation of two states can be seen as a process in which 
the first party (Alice) encodes a direction into a quantum state while the receiving party 
(Bob) attempts to estimate the signal without having a classical reference frame relative to 
which he can measure it. Therefore Bob resorts to finding the relative orientation of the 
signal state with respect to the orientation of some given state (say, a coherent state of spin 
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j 2 ), which serves as a quantum reference frame. So far the value of j 2 has been kept equal 
to 1/2. We move on to consider what happens as we increase the value of j 2 . The limit 
j2 - * oo could be regarded as the limit in which the quantum reference direction becomes a 
classical one. 

As before, we need to consider states with mi = and mi = 1. Different blocks of p(f3) 
are found by coupling the representations of \ipi) and \i/j 2 )- For the mi = 1 sector, i.e. the 
case of parallel spins, p(/3) has three one- dimensional blocks, since 



The optimal measurement is then given by projections on subspaces of total angular mo- 
mentum J. 

For the m x = sector we have 



thus p(f3) and the operators have two 1-dimensional blocks and one 2-dimensional block. 
The POVM elements acting on the J = j 2 — 1 and J = j 2 + 1 subspaces are rank one 
projectors onto these subspaces. In the J = j 2 subspace, an optimization procedure similar 
to the one in Sec. IIHAI needs to be done. 

Carrying out the optimization, we find that the optimal state | ipopt) belongs to the mi = 
sector for all values of j 2 . The optimal state is not fixed but rather depends on the value of 
j 2 . The limit j 2 — > oo yields the optimal asymptotic state 



with doo = 0.595. The dependence of the state | ip op t) on j 2 is plotted in Fig|21 

Comparing the fidelity achieved with the optimal state, as a function of j 2 , to the fidelity 
obtained using both the anti-parallel and the parallel spin states, we get, quite remarkably, 
almost the same fidelity for any value of j 2 . The parallel spins give a slightly higher fidelity, 
with the difference (already very small for j 2 = 1/2) rapidly decreasing for increasing values 
of j 2 . This comparison is plotted in Fig which shows that the plot for the parallel and 
anti parallel spin states coincide. 

Notice that although we expect the limit j 2 — ► oo to be equivalent to measuring the 
state | ^i) against a classical reference direction, the anti-parallel spin states do not become 



1 ® h = {32 - 1) © J2 © C?2 + 1). 



(0 © 1) © j 2 = j 2 © (j 2 - 1) © j 2 © 0" 2 + 1), 




(62) 
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FIG. 2: The optimal state parameter a (vertical axis) with increasing values of j2 (horizontal axis). 

Fidelity 




Log j 



FIG. 3: Fidelity achieved using the optimal state (red, upper line), anti-parallel spins and parallel 
spins (green, lower line) as a function of jV 

optimal in this limit. This seems to be in contradiction with the result of 0, |fj], who 
showed that if only two spins are available, the optimal direction indicator is provided by 
anti-parallel spins along that direction. The resolution to this apparent contradiction relies 
on the fact that the jo — > oo limit considered here is not equivalent to an estimation of a 
direction with respect to a classical reference frame, as , but rather to 

that of an angle between the same vector, and, say, the z-axis of such a frame. In the next 
section we will consider the latter estimation task and show that it coincides with the limit 
j2 oo discussed here, demonstrating that a macroscopic spin can be treated as a classical 
reference direction. 

It is conceivable to consider a different estimation task in which one observer, say Alice, 
indicates a direction in space using the state while Bob encodes his reference frame 

into the state \ip2), and finally a third observer is interested in the orientation of Alice's 
direction in Bob's frame. This task, however, cannot be performed using a spin coherent 
state | ^2) = |j2, J2) ■ It would be interesting to see what state would encode Bob's frame 
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in an optimal manner, and whether taking the appropriate limit would reproduce previous 
results for direction alignment. 

C. Quantum-classical correspondence 

Let us consider a scenario in which Alice prepares a state, in order to indicate a 

chosen direction, and Bob is requested to estimate the angle between this direction and the 
z-axis of his classical reference frame, which replaces the quantum reference direction (| 
in the previous study. We assume no knowledge of Alice's reference frame, and without loss 
of generality we can assume that Alice chooses to indicate her z-axis. If the transformation 
relating Alice's frame to Bob's is parameterized by the Euler angles x> P an d 4>, then, in the 
latter reference frame, Alice's state is given by U(x, P, <p)\ipi). Bob's task is now to estimate 
the angle ft between Alice's z-axis and his z-axis. As before, given l^i), we are seeking a 
POVM that maximizes the fidelity 

F{E,} = Y. I d x ^Tr[a 1 ( X ,PA)E,}f(^P), (63) 

where d x ^ = sin ftdpd<j)dx/87r 2 is the invariant measure of the rotation group and 

a x ( X , P, 0) = U{ X , P, 0) (i>iMx, P, 

Since we are not interested in estimating the angles <ft an d X> we can integrate over them 
and define an averaged density matrix function of P only, 

°(P) = ^ j d(f>dxai{x,PA)- (64) 

Let us analyze the form of the density matrix a(P). Writing | ipi) as in Eq. (|14|) . and using 
the definition of the rotation operator matrix elements 

(jm'\U( X ,PA)\jm) = e im '*dl, m (PY m t (65) 

we derive the matrix elements of <j(P) as 

(j'm'\a(P)\jm) = ^5 m / m a^'af d?^, r (P)d° mr {P) ■ (66) 

r 

The matrix a(P) is thus diagonal in the indices m, but has off-diagonal elements with 
different values of j. In this respect it is similar to the matrix p(j3) defined in Eq. (j3*Uj) . 
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which was diagonal in the representation J with off-diagonal elements between different 
values of j±. Indeed, by taking the limit 22 00 in Eq. (JHUj) and interchanging the indices 
J — > m, we get the asymptotic equivalence between the matrices p{0) and a((3), i.e., 

lim (JjMPWi) = (jim\a(P)\f im ). (67) 

This result shows that the fidelity achieved by relating any state | ^1) to a classical reference 
direction is identical to the one achieved in the limit j'2 —* 00 discussed in the previous 
section (with the same state \ipi))- Consequently, also the optimal state will be identical in 
the two cases. 



IV. CONCLUDING REMARKS 



In this work we studied the problem of estimating relative rotation angles of quantum 
signals. By considerations of global symmetry invariance, we derived the general form of 
the signal state and of the appropriate set of measurements and estimation strategies. For 
special low-dimensional cases, explicit optimization is carried out. 

With these tools we have studied the state preparations which maximize the average 
fidelity of the estimation procedure, and compared these with the results found in the esti- 
mation of absolute rotations. We have also discussed the asymptotic limit in which one of 
the quantum states becomes a macroscopic spin. In this limit, the resulting estimation task 
is identical to an estimation of the orientation of a quantum state with respect to a classical 
reference direction. 

Many important questions remain open for investigation. A broader extension of the 
problem would lead to the estimation of the orientation of a quantum state relative to 
another one that encodes a full reference frame (three axes), rather then a single direction. 
In such a scenario one has to estimate two angles (polar and azimutal). An even more 
elaborate framework is one in which each quantum state encodes a full reference frame, 
and one is interested in estimating the transformation that would align these two reference 
frames. It would be interesting to compare the optimal states found in all these cases with 
those found in prior studies Iflflfll, when one of the reference frames i 8 classical. We 
would like, finally, to emphasize that we did not include in our communication scheme the 
possibility of the two parties sharing a prior entangled state. In this case, the encoding and 
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the detection of relative information could proceed via a covariant dense coding scheme, 
which could increase the efficiency of the estimation. 
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